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ABSTRACT 

We have calculated the rate coefficients for D(ls) + H + ^ D + + H(ls) using 
recently published theoretical cross sections. We present results for temperatures 
T from 1 K to 2 x 10 5 K and provide fits to our data for use in plasma modeling. 
Our calculations are in good agreement with previously published rate coefficients 
for 25 < T < 300 K, which covers most of the limited range for which those 
results were given. Our new rate coefficients for T > 100 K are significantly 
larger than the values most commonly used for modeling the chemistry of the 
early universe and of molecular clouds. This may have important implications 
for the predicted HD abundance in these environments. Using our results, we 
have modeled the ionization balance in high redshift QSO absorbers. We find 
that the new rate coefficients decrease the inferred D/H ratio by < 0.4%. This is 
a factor of > 25 smaller than the current > 10% uncertainties in QSO absorber 
D/H measurements. 

Subject headings: atomic data - atomic processes - early universe - interstellar: 
molecules - quasars: absorption lines 



1. Introduction 

Deuterium plays an important role in addressing several fundamental questions in astro- 
physics. The deuterium abundance is a key constraint for models of big bang nucleosynthe- 
sis. Primordial D /H measurements provide the most sensitive probe of the baryon-to-photon 
density ratio rj. This, in combination with the cosmic microwave background measurement 
of the photon density, can be used to determine the cosmological baryon density (Buries & 
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Tytler 1998; Lemoine et al. 1999; Tytler et al. 2000). Deuterium may also be important 
in the formation of structure in the post-recombination era of the early universe. HD is 
the second most abundant primordial molecule, after H 2 , and cooling radiation from it may 
play a role in the formation of the first collapsing objects (Puy et al. 1993; Stancil, Lepp, & 
Dalgarno 1998). Lastly, as the universe evolves, deuterium is slowly destroyed in stars where 
it is burned into 3 He. Mapping the temporal and spatial variations in the D/H abundance 
ratio can shed light on the time history of star formation in different regions of the cosmos 
(Tosi 1998; Tosi et al. 1998; Lemoine et al. 1999). 

Investigations into these issues are carried out through studies of gas phase D and 
deuterated molecules in, for example, the early universe (Galli & Palla 1998; Stancil et 
al. 1998), QSO absorption systems (O'Meara et al. 2001), molecular clouds (Tielens 1992; 
Bertoldi et al. 1999; Wright et al. 1999) and the interstellar medium (Linsky et al. 1995). 
Stancil et al. (1998) have given a recent listing of deuterium models for the early universe, 
molecular clouds, and the interstellar medium. Recent reviews of deuterium observations 
have been given by Lemoine et al. (1999) and Tytler et al. (2000). 

Interpreting these studies requires an accurate knowledge of all collision processes in- 
volving D. Particularly significant is the near-resonant charge transfer (CT) process 

D(Ls) + H+ -> D+ + H(ls) (1) 

and the reverse process 

D+ + H(Ls) -> D(Ls) + H+ (2) 

In the early universe these are two of the most important processes involving deuterium 
(Galli & Palla 1998). In molecular clouds, Process 1 followed by the exothermic reaction 

H 2 + D+ -> HD + H+ (3) 

is a major source of HD (Black & Dalgarno 1973; Dalgarno, Black, & Weisheit 1973; Watson 
1973). This is also an important source of HD in the early universe (Stancil et al. 1998). 

Recently Igarashi & Lin (1999) and Zhao, Igarashi, & Lin (2000) have carried out 
cross section calculations for Reaction 1 (o^) and Reaction 2 (<r 2 ) using a hyperspherical 
close-coupling method. This technique is free from the ambiguities associated with the 
conventional Born-Oppenheimer approach. Here we use their results to produce CT rate 
coefficients for Reactions 1 (oji) and 2 (0:2)- In § 2 we describe how we evaluated ai and 0:2- 
Our results are presented in § 3 and compared with previously published calculations. Some 
astrophysical implications are discussed in § 4. 
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2. Calculation of the Rate Coefficients 

We consider only capture from and into the Is level of H and D. The reason for this is 
twofold. First, in the sources discussed in § 1, neutral H and D are expected to be found 
essentially only in their ground state. Second, at the low temperatures (i.e., low collision 
energies) relevant for these sources (T < 3 x 10 4 K [i.e., k B T < 3 eV]), CT into the Is level 
is predicted to be over 4 orders of magnitude greater than capture into other levels. CT 
into higher levels does not become important until collision energies of > 10 3 eV (Alt et al. 
1994). 

We use the results of Igarashi & Lin (1999) and Zhao et al. (2000) for ai(E) and a 2 (E) at 
center-of-mass energies E from 2.721 x 10~ 8 to 2.721 eV (Igarashi 2001). Due to the binding 
energy of D(ls) being slighlty larger than that of H(ls), Reaction 1 is endothermic with a 
threshold of 3.7 meV (43 K). Hence, o\ is predicted to be smaller than a 2 at all energies, but 
particularly at low energies. As E increases, <j\ and a 2 converge, and for E > 2.72 eV, the 
two are predicted to lie within < 0.1% of one another (Zhao et al. 2000; Igarashi 2001). The 
energy dependences for o\ and o 2 allow us to extend the results of Igarashi and collaborators 
to higher energies. We do this using the calculated cross sections of Dalgarno & Yadav 
(1953) for the related reaction 

H+ + H(ls) -> H(ls) + H+. (4) 

First, to extend the data for a 2 , we multiply the energy scale of Dalgarno & Yadav by 
Hhd/^hh where \i is the reduced mass for the HD and HH systems. This effectively matches 
the velocity scale for each data set. Then, we multiply the results of Dalgarno & Yadav by a 
factor of 0.959 to set it equal to the results of Igarashi and collaborators at 1.333 eV. Next, 
for energies between 1 and 2.72 eV, we fit the ratio of 0\jo 2 from Igarashi to the formula 

which yields A = 2.692 x 10~ 3 and B = 7.936 x 10~ 4 . We calculate o\ for energies above 
« 2.72 eV using the scaled cross sections of Dalgarno & Yadav multiplied by Equation 5. 

We use the resulting data for o\ and a 2 from energies of 2.721 x 10~ 8 to 10 3 eV to 
evaluate the rate coefficients a±(T) and a 2 (T) as a function of the gas temperature T. 
Rate coefficients are calculated numerically using the desired cross section times the relative 
velocity and convolving these results with the appropriate Maxwellian distribution (taking 
the reduced mass into account). Cross sections for energies not calculated by Igarashi and 
collaborators or by Dalgarno & Yadav are evaluated using a spline interpolation method 
(Press et al. 1992) for a(E) versus \og(E). 
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3. Results and Comparisons 

Our calculated results for a\(T) and 0:2 (7") are given in Tab. 1 from 1 K up to 2 x 10 5 K. 
These results are also plotted in Fig. 1 from 1 to 30,000 K. As expected the rate coefficient 
for the endothermic Reaction 1 decreases dramatically for T < 10 K. The rate coefficient for 
Reaction 2 decreases slowly with decreasing temperature down to T pa 15 K. Below this a 2 
begins to increase with decreasing temperature. We attribute this to the rapid increase in 
<7 2 with decreasing collision energy (see Fig. 2 of Igarashi & Lin 1999). 

We have fitted our calculated CT rate coefficients using the formula 

a(T) = aT b exp(-c/T) + dT e . (6) 

The best fit values are listed in Tab. 2. The fits for ot\ and a 2 are accurate to better than 
6% and 4%, respectively, for 2.7 K < T < 2 x 10 5 K. 

Several other groups have carried out detailed calculations for o\ and a 2 . Davis & 
Thorson (1978) published results for o\ from 3.7 to pa 100 meV. Results for o\ and o 2 
were reported by Hunter & Kuriyan (1977) for energies from 10~ 3 to 7.5 eV, by Hodges 
& Breig (1993) from 10~ 3 to 10 eV, and by Esry et al. (2000) from pa 3.7 meV to 8 eV. 
In general the calculations of Igarashi and collaborators are in good to excellent agreement 
with these published results. The most significant difference is for energies above 10" 3 eV 
where the results of Hunter & Kuriyan can fall as much as 15% below those of Igarashi and 
collaborators (Igarashi & Lin 1999; Zhao et al. 2000; Igarashi 2001). This is partially due 
to the accidental overlap of minima in the oscillating cross sections with the energy points 
published by Hunter & Kuriyan (1977). 

There have been a couple of experimental measurements of o\. We are unaware of 
any experimental results for a 2 . Absolute measurements of (j\ have been carried out by 
Newman et al. (1982) for energies between pa 0.1 and 10 eV. The theoretical results of 
Zhao et al. (2000) and Esry et al. (2000) are in good agreement with these measurements. 
Relative measurements for o\ have carried out by Wells et al. (2001) for energies between 
threshold and 1 eV. Good agreement was found with the calculations of Esry et al. (2000) 
between pa 0.02 and 1 eV. Uncertainties in background subtraction limit the reliability of 
the experimental results below 0.02 eV. 

Using the results of Hunter & Kuriyan (1977), Watson et al. (1978) calculated a± and 
a 2 for a number of temperatures between 10 and 300 K. These results are listed in Tab. 1 
and also plotted in Fig. 1. For T > 50 K, the results of Watson et al. agree with ours to 
better than 5%. At 25 K their results differ from ours by pa 8% and at 10 K by pa 30%. The 
differences for T < 25 K are most likely due to the uncertainty associated with extrapolating 
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the results of Hunter & Kuriyan to energies below those published (Watson et al. 1978). 

Galli & Palla (1998) fit the results of Watson et al. (1978) for a,\ and a 2 . The resulting 
fitted rate coefficients are plotted in Fig. 1. Between 10 and 300 K, these fitted rate coef- 
ficients agree with our results not quite as well as the results of Watson et al. Agreement 
with our results becomes progessively worse the further one extrapolates these fitted rates 
outside this temperature range. 

Watson (1976) presented an estimate for cci and a 2 which we plot in Fig. 1. These 
estimated rate coefficients are in poor agreement with our results here, differing significantly 
in both the values and temperature dependences of a± and a 2 . 

4. Some Astrophysical Implications 

4.1. The Early Universe 

Recently, results from a number of different chemical models of the early universe have 
been published. For these models, Puy et al. (1993) and Standi et al. (1998) used the 
estimated rate coefficients of Watson (1976). Galli & Palla (1998) used their fits to the 
results of Watson et al. (1978). For redshifts z > 50, where the gas temperature is predicted 
to be > 50 K (Puy et al. 1993), the rate coefficients used by Puy et al. and Stancil et al. 
begin to differ significantly from our newly calculated results. At z 400 (T 1000 K; Puy 
et al. 1993), the rate coefficients used by them are a factor of rs 3 smaller than our results. 
In contrast, the extrapolated rate coefficients used by Galli & Palla are only xs 22% smaller. 
Determining the full implications of our new rate coefficients will require re-running updated 
versions of these various chemical models of the early universe. 

4.2. Molecular Clouds 

Modeling studies of molecular clouds have been carried out recently by Millar, Bennett, 
& Herbst (1989), Pineau des Forets, Roueff, & Flower (1989), Heiles, McCullough, & Glass- 
gold (1993), Rodgers & Millar (1996), and Timmermann (1996). These studies have all used 
the results of Watson (1976) for a± and a 2 and hence significantly underestimate these two 
rate coefficients for T > 100 K. Because Reaction 1 followed by Reaction 3 is predicted to 
be a major source of HD in molecular clouds (Black & Dalgarno 1973; Dalgarno, Black, & 
Weisheit 1973; Watson 1973), underestimating a\ could in turn lead to an underestimate in 
the amount of HD produced in these clouds. 
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4.3. High Redshift QSO Absorption Systems 

Observations of high redshift QSO absorption systems are used to infer the primordial 
D/H ratio. These studies are carried out assuming that the D I/H I ratio is identical to that 
of D/H (Buries & Tytler 1998). Here we investigate the validity of this assumption in light 
of the different values for ot\ and a 2 - 

The D/H ratio inferred from these observations is given by 

n D f Ul N(Di) 



n H /diN(Hi) 



(7) 



where n D is number density of D, N(D i) is the column density of D I, and / Dl is the 
abundance of D I relative to the total abundance of D. Similar definitions exist for H and 
H I. We can write / D t as 

^- = l + =2=. (8) 

/Di "Dl 

A similar expression can be written for /hi- 

Currently there are believed to be six reliable measurements of D/H in high redshift 
QSO absorbers (Pettini & Bowen 2001). These measurements all assume /bi = /hi- To 
determine the validity of this assumption, we evaluate 

/hi _ 1 + »dii/wdi ^ 
/ D i 1 + n H n/ i 

using our new results for a.\ and a 2 - 

At the inferred temperatures in these six absorbers (T ps 1.1 x 10 4 K), the gas is predicted 
to have an insignificant abundance of molecules (Petitjean, Srianand, & Ledoux 2001). The 
ionization balance of D in these QSO absorbers can therefore be written 

ii _ Ad i + n c C Dl + nHiiQp + H+ + n xi+®p+xi+ 
n D i n c R D „ + n H ia H +D+ + n x«x+D+ 

In the numerator on the right-hand-side (RHS) of this equation, /3d i is the photoionization 
(PI) rate of D I due to the radiation field and accounts for further ionization due to the 
resulting non-thermal photoelectrons, n e is the electron density, Cd l is the electron impact 
ionization (EII) rate coefficient due to thermal electrons, a D+H + is the rate coefficient for 
D + H + collisions producing D + , and similarly for a D+X 9+ where X q+ represents a g-times 
charged ion of element X and the sum over X q+ includes the ions of all elements except 
for those of H and D. Here, q;d+h+ and a D+ x«+ are purely CT rate coefficients. This is 
because collisions which leave both colliding particles in an ionized state are predicted to be 
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insignificant at the temperatures of interest (Janev, Presnyakov, & Shevelko 1985). In the 
denominator on the RHS, -Ron is the radiative recombination (RR) rate coefficient for D II, 
Q! h+d+ is the CT rate coefficient for H + D + collisions producing D, and similarly for «x+d+- 

The PI rates and EII and RR rate coefficients for D and H are expected to be essentially 
identical (Galli & Palla 1998; Standi et al. 1998). The differences in the energy level structure 
of D and H have an insignificant effect on these processes. Thus in Eq. 10, we can substitute 

Phi = Pdi, (11) 

Chi = Cdi, (12) 

and 

-Rhu = -Rdh- (13) 

Next we add and subtract 

KB n«H+D+ + ^2 n X9+«H+X9+ (14) 

to the numerator on the RHS of Eq. 10 and 

riB i«d+h+ + ^2 n x«x+H+ (15) 

X 

to the denominator. We note that using Eq. 10 we get /inn/Vim by interchanging all charge 
states of D with the corresponding charge states of H (and vice versa). Hence, we can rewrite 
Eq. 10 as 

»D ii _ «Hn ( 1 ±Jl/8 \ 
"Di "Hi \1 +I2/8J 



where 



(16) 

7l = « D+ H+ "H+D+ + > , "D+X9+ - «H+X-?+ , (17) 

WHii ^ WHii V / 



72 = «H+D+ — «D+H+ + — f «X+D+ - "X+H+ J , (18) 

nu i ^ n H i V / 



and 



= Kru H "D+H+ + > "X+H+- (19) 



%i "Hi X UH 1 



We can simplify 7x and 72 . First we note that «d+h+ = a i and q;h+d+ = ®2] and at 
the temperatures of interest a,\ ~ «2 ~ 8.3 x 10 -9 cm 3 s" 1 . Now, to a first approximation 
ft-Dii/ft-Hn and wdi/ihi will be equal to the primordial D/H value which we take to be 
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2 x 10 5 (from Pettini & Bowen 2001). Hence the second term in Eqs. 17 and 18 is 
roughly 5 orders of magnitude smaller than the first term and can be dropped. 

At energies important for 1.1 x 10 4 K, we note that U\ a 2 . Similarly, we expect 
at these temperatures <7b+xi+(v) ~ cr H+X g+(i>) and <tx+d+(v) ~ o"x+h+(^)) where v is the 
relative velocity. As a result, we have a D+ xi+ ^Hdx ~ <2h+x«+ y/^HX and a x +D+ y/J^DX ~ 
«x+h+ \[Whx- Here, ~ 1 and /i£>x ~ 2 are the reduced masses. For those ions where CT 
is important in photoionized plasmas (e.g., QSO absorbers) we estimate that aD+x<?+, o;h+X9+, 
«x+d+, and «x+h+ will all be < 10~ 9 cm 3 s _1 (Kingdon & Ferland 1996). The expressions 
in the parenthesis in the third term in Eqs. 17 and 18 are thus < 4 x 10~ 10 cm 3 s" 1 . 
Furthermore, we note that the metallicity in these absorbing systems are ~ 10~ 2 solar 
(Pettini & Bowen 2001) and that we expect %+/nHn and %/nHi will be within a couple 
of orders of magnitude of these reduced abundances. As a result, we can also drop the third 
term in Eqs. 17 and 18. With these approximations we can rewrite Eq. 16 as 

n D n n Uu f aA fl + 5/aA 

111 1 (20) 



As «i and a 2 differ by < 0.4% at T « 1.1 x 10 4 K, Eq. 20 reduces to 

^I^W. (21) 

«Di "-Hi V a 2/ 

Because <Ji has a threshold, in contrast with <t 2 , the ratio ai/a 2 will always be less than 1. 
For the temperature of interest here, we have 

I^Eli ~ 0.996^. (22) 
n D i n H i 

Substituting this results into Eqs. 7 and 9, we find 

n D / « 0.996A^(Di)/A^(Hi) (n H „>n Hl 



n H I 1.000iV(D i)/JV(H i) (n H „<n Hl ). ^ 

The uncertainty in the above factor of 0.996 depends on the accuracy of the theoretical cross 
sections we have used here. Verifying the accuracy of these cross sections will require further 
theoretical and experimental studies. However, we note that the variation of m 0.4% in 
nc/«H is a factor of > 25 smaller than the current > 10% uncertainties in QSO absorber 
D/H measurements. Thus, it is likely to be some time before the subtle differences in 
Reactions 1 and 2 become important for primordial D/H measurements. 
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Fig. 1. — Calculated rate coefficients for cci(T) [D(ls) + H + — > D + + H(ls)] and the reverse 
process a 2 (T) [D + + H(ls) — > D(ls) + H + ] versus gas temperature T. Results for a.\ are 
from this work {thick solid curve), Watson (1976; thick dashed curve), Watson et al. (1978; 
squares), and Galli & Palla (1998; thick dotted curve). Results for a 2 are from this work 
(thin solid curve), Watson (1976; thin dashed curve), Watson et al. (1978; circles), and Galli 
& Palla (1998; thin dotted curve). 
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Table 1. Calculated rate coefficients for «i(T) [D(ls) + H + — > D + + H(ls)] and the 
reverse process a 2 (T) [D + + H(ls) — > D(ls) + H + ] versus gas temperature T. 
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Table 1 — Continued 
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8.37E-09 


12000 


8.62E-09 


8.65E-09 


13000 


8.89E-09 


8.92E-09 


14000 


9.15E-09 


9.17E-09 


15000 


9.39E-09 


9.42E-09 


oneiric* 
ZUUUU 


l.Uoriv-Uo 


l.Uoriv-Uo 


25000 


1.14E-08 


1.14E-08 


30000 


1.22E-08 


1.23E-08 


35000 


1.30E-08 


1.30E-08 


40000 


1.37E-08 


1.37E-08 


50000 


1.49E-08 


1.49E-08 


75000 


1.73E-08 


1.74E-08 


100000 


1.93E-08 


1.93E-08 


200000 


2.50E-08 


2.50E-08 



a Watson et al. (1978) 



Table 2. Fit parameters for our calculated rate coefficients Oi\(T) 
[D(ls) + H+ -> D+ + H(ls)] and the reverse process a 2 {T) [D+ + H(ls) -»• D(ls) + H+]. 



Rate Coefficient 


a 


6 


c 


d 


e 




(cm 3 s _1 ) 




(K) 


(cm 3 s _1 ) 






2.00E-10 


0.402 


37.1 


-3.31E-17 


1.48 


« 2 


2.06E-10 


0.396 


33.0 


2.03E-09 


-0.332 



